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Abstract 

In this paper, we consider probability measures and v ona d-dimensional 
sphere in R'^"'"^, d > 1, and cost functions of the form c(x, y) = ^(^^-j-^) that 
generahze those arising in geometric optics where l{t) = — logt. We prove that 
if ^ and u vanish on {d — l)-rectifiable sets, if \l'{t)\ > 0, hm(_^o+ l{t) = +00, 
and g{t) := t{2 — t){l'{t))'^ is monotone then there exists a unique optimal 
map To that transports /i onto i^, where optimality is measured against c. 
Furthermore, infx jToX — x| > 0. Our approach is based on direct variational 
arguments. In the special case when l{t) = — logt, existence of optimal maps 
on the sphere was obtained earlier in [8] and [22] under more restrictive as- 
sumptions. In these studies, it was assumed that either ^ and v are abso- 
lutely continuous with respect to the ^-dimensional Haussdorff measure, or 
they have disjoint supports. Another aspect of interest in this work is that 
it is in contrast with the work in [7j where it is proved that when l{t) = t 
then existence of an optimal map fails when /x and v are supported by Jordan 
surfaces. 
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1 Introduction 

In Euclidean space R"^"*"^ consider a reflector system consisting of a point source O 
radiating with intensity /(x) in directions x G X, where X is a closed set on a 
dimensional unit sphere S*^ ^ R'^"'"^ centered at O, and a smooth perfectly reflecting 
hypersurface R, star-shaped relative to O, which intercepts and reflects the light rays 
with directions from X; see Fig. [H Assuming the geometric optics approximation 
and applying the classical reflection law to determine the set of reflected directions 
Y C S'^ (after one reflection), we obtain an associated with R "reflector map" 
^ : X — > Y. Assuming that no energy is lost in the process, one can apply the 
energy conservation law to calculate the intensity distribution L{y) produced on Y. 
The reflector problem consists in solving the inverse problem in which the source 
O, the sets X, Y and the intensities / and L are given in advance and the reflector 
R needs to be determined. That is, R should be such that ^(X) ^ Y and 

^(e(x))|J(0(x)|=/(x) 

for all X in the interior of X; here J(^) denotes the Jacobian determinant of ^. 

Problems of this type arise often in applications, for example, in design of re- 
flector antennas [23]. In various forms the reflector problem has been considered 
by many authors and numerous papers by engineers (at least since early 60-th [12] 
until now [22]) are devoted to this subject. It was introduced in electrical engineer- 
ing and optics independently of the mass transport problem. Because of the strong 
nonlinear constraints appearing in the problem, progress has been slow and many 
theoretical and computational issues still remain open. The problem continues to 
attract considerable attention because of its practical importance and mathematical 
subtleties. It may be pointed out that a version of the reflector problem appears on 
the famous list of unsolved problems proposed by S.T. Yau [21] in 1993. 

Analytically, the reflector problem considered in this paper can be formulated as 
a nonlinear second order elliptic partial differential equation of Monge- Ampere type 
on a subset of S'^. In such form it has been studied by V. Oliker and P. Waltman 
[16], L. Caffarelh and V. Oliker [3], X.-J. Wang [2T], P.G. Guan and X.-J. Wang 
|lUj . L. Caffarelli, S. Kochengin and V. Oliker [4], V. Oliker [TTj, and other authors. 

Recently, T. Glimm and V. Oliker |8j and, independently, X.-J. Wang [22j have 
shown that if the function I (resp. L) is treated as the density of a measure /i (resp. 
z/) that are absolutely continuous with respect to the volume measure on S'*, then 
the reflector problem can be studied as a variational problem in the framework of 
Monge-Kantorovich theory, that is, a problem of finding an optimal map minimizing 
the transport cost of transferring fi onto u with the cost function — log(l — x ■ y). 
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In contrast with other cost functions considered usually in the Monge-Kantorovich 
theory, this cost function may assume infinite values. Consequently, in order to 
overcome this difficulty, a geometric condition requiring the supports of fi and u 
to be disjoint was imposed in [8] and [22] to establish existence and uniqueness of 
optimal maps. Without imposing the condition that the supports of the measures 
are disjoint, existence and uniqueness of optimal maps was also obtained in [S]. 
However, the proof is indirect as it relies on existence of weak solutions in the 
reflector problem established earlier in [3], [17]. 

The contribution of this study is twofold. First of all, we obtain existence and 
uniqueness of optimal maps To for a class of cost functions that may be infinite. This 
class includes the logarithmic cost function of the reflector problem as a special case. 
The cost functions are precisely of the form c(x, y) = ) where \l'{t)\ > 0, 

limt^o+ K^) = +00, and g(t) := t{2 — t){l'(t)Y is monotone. Furthermore, we prove 
that infx |ToX — x| > 0, which we view as an intermediary step in the study of 
the regularity of the map Tq. Secondly, our approach is variational and direct; the 
supports of measures /i and u are allowed to overlap and it is merely required that 
these measures vanish on {d — l)-rectifiable subsets. The precise statement can be 
found in Theorem 14. 6[ 

Let us recall the main principles which ensure existence of optimal maps, which 
apparently, have been explicitly pointed out for the first time in [5] and later ex- 
ploited by many authors. Assume we are given a cost function c : R'^+^ x R'^+i ^ R 
and two probability measures fi and u on R'^+^, say, absolutely continuous with 
respect to Lebesgue measure. Existence and uniqueness of an optimal map trans- 
porting /i onto u against c is ensured if Vxc(x, •) is a one-to-one map of R'^"'"^ into 
itself. Note that if c(x, y) = |x — yp then Vxc(x, •) has this property. 

It is shown in [7] that if instead, /x and z/ are supported by S*^, then existence of 
an optimal map transporting /i onto u against c(x, y) = |x — yp, fails. One of the 
points in our study here is that we start with a cost function c(x, y) = Z(|x — yp/2) 
such that Vxc(x, ■) fails to be a one-to-one map of R*^"*"^ into itself. However, if 
x G S*^ and a G R^^^ is distinct from the origin, we observed that there exists 
a unique y G satisfying the equation Vg<jc(x, y) = a. Here, VgdC(x, ■) is the 
tangential gradient of c on S'* with respect to x. The injectivity of VgdC(x, ■) yields 
existence of an optimal map that transports onto u, if /x and z/ are supported by 
S*^ and vanish on (rf — l)-rectifiable sets. We use the Kantorovich theory as a tool to 
give a transparent explanation to the existence of a unique solution in the reflector 
problem, under sharp assumptions. Our main results are stated in Theorem 14.61 
We refer the reader to a variant of the reflector problem involving two reflectors 
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considered in a paper by T. Glimm and V. Oliker [9]. We also refer the reader to a 
recent study by N. Ahmad [2] in the plane, motivated by [7]. 

The paper is essentially self-contained. It is organized as follows. In order to 
motivate out subsequent considerations, we begin with a review of the reflector 
problem in section [2l In section [3] we review and extend some results from the 
Monge-Kantorovich theory. Our main results establishing existence and uniqueness 
of optimal maps are in section |H 

2 A review of the reflector problem 

Let X, Y, /, L and R be as in the introduction. If n is the unit normal field on 




Figure 1: The reflector problem 

i?, then the incident direction x and the reflected direction y are related by the 
reflection law 

y = x — 2(x-n)n. (1) 

Thus, the hypersurface R defines the reflector map ^ : x ^ y which maps the "input 
aperture" X C S'^ onto the "output aperture" Y C S'^; see Fig. 1. The intensity of 
the light reflected in direction y = ^(x) is given by /(x)/| J(^(x))|. 

Suppose now that the closed sets X and Y on S'^ are given as well as nonnegative 
integrable functions / on X and L on Y representing, respectively, the intensity of 
the source and the desired intensity on the far-region Y. The reflector problem is 
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to determine a reflector R such that the map ^ defined by R maps X onto Y and 
satisfies the equation 

L{y) = l{r\y))\J{r\y))l y g MY; (2) 

see [15], [16]. Note that this is an equation on the output aperture Y rather than on 
the input aperture X. One could also set it up on X [14j , but ([2]) is more convenient 
for our purposes here. 

It was shown in [151 US] that there exists a scalar quasi-potential p : Y —>■ [0, oo) 
from which the reflector R can be recovered and in terms of which the equation ([2]) 
when J{^^^) 7^ is the following equation of Monge- Ampere type 

L(y) = /(rHy)) '''''^";t^+';"^'^'' . ^'^■"'Y- 

where e is the standard metric on S'^, p = {p"^ + |Vpp)/2p, and Hess{p), Vp are 
computed in the metric e. In terms of p the position vector of R is given by 

r(y) = -Vp(y) - (p(y) - p(y))y : Y R^+\ (4) 

while 

^-1/ ,x _ r(y) 



P(y) 

Note that |r| = p. 

A close examination of (jl]) shows that it describes R as an envelope of a family 
of paraboloids of revolution P{y) tangent to R, parametrized by their axes y G Y 
and with common focus O. For each y, p{y) is the focal parameter of -P(y). This 
observation was used in [3] for the weak formulation of the reflector problem where 
a class of convex reflectors corresponding to positively elliptic solutions of ([3]) was 
introduced. Reflectors corresponding to negatively elliptic solutions of can be 
introduced and analyzed by similar methods [S], [22] • Such reflectors, however, are 
only piecewise concave (relative to the origin O). For brevity we discuss here only 
convex reflectors which we now define. 

Let Y be a subset on S*^ and p' : Y — > (0, oo) a bounded function. Let {-P(y)} 
be a family of paraboloids of revolution with axes of direction y G Y, common focus 
O and polar radii 

Py(x) = Y^^, X G S'^ \ {y}. (5) 
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The closed convex hypersurface R given by r(x) = p(x)x, x G S**, with 

p(x) = inf Pj;(x), X G S*^, (6) 

yeY 

is called a reflector deflned by the family {-P(y)}yGY (with the light source at 
O). 

Let i? be a reflector, z a point on i?, P a paraboloid of revolution with focus at 
O and B the convex body bounded by P. If P C P and z G P then P is called 
supporting to R at z. 

For any y G S'^ a reflector R has a supporting paraboloid with axis y. The 
corresponding continuous function giving the focal parameters of all such supporting 
paraboloids P(y), y G S*^, of R is called the focal function and denoted by p. The 
natural question of characterization of focal functions of closed convex reflectors was 
partially answered in [17j. 

The reflector map (possibly multivalued) generalizing ([1]) and denoted again 
by ^ is deflned for x G S*^ as 

^(x) = {y G S*^ I P(y) is supporting to R at p(x)x}. (7) 

It follows from ([5]) and ([6]) that, equivalently, the reflector map can be deflned as 

e(x) = {y G I logp(x) - logp(y) = - log(l - x • y)}. (8) 

The total amount of energy transferred from the source O ma, given set of directions 
cij C Y is best described with the use of the map deflned on any subset C Y 
by setting 

It is shown in [3[ [T7] that for any Borel subset uj G Y the set is Lebesgue 

measurable on S*^. Thus, if / > is the intensity of the source then the total amount 
of energy transferred by R to the set u; C Y is given by the "energy" function 

G{R,u;) = [ /(x)da(x), (9) 

where da is the standard rf— volume form on S'^. It is assumed here and everywhere 
below that I is extended from X to the entire by setting /(x) = Vx G 
S*^ \ X. The function G{R, u) is a Borel measure on Y (not necessarily absolutely 
continuous) [5], [T7]. 
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For a given nonnegative and integrable function L on Y we say that a closed 
convex reflector i? is a weak solution of the reflector problem if 

G{R,uj)= L{y)da{y) for any Borel set uCY. (10) 

Put 



^\A\ = / /(x)da(x), u[B] = / L{y)daiy) (11) 

J A J B 

for y4 C X and i? C Y Borel sets. An obvious necessary condition for existence of 
a weak solution to the reflector problem is that the total energy of the source and 
the total energy on the output aperture are in balance, that is, 

/x[X] = y[Y]. (12) 

Excluding the trivial case when either /i[X] or z/[Y] is zero, it may be assumed, 
without loss of generality, that measures satisfying ( fT2l) are normalized so that they 
are probability measures. 

The following existence result was established by L. Caffarelli and V. Oliker in 
^ (and reproduced partly in [1] and [IZ]). 

Theorem 2.1. Let X and Y he closed sets on S** (possibly coinciding with S'^) and 
I > 0, L > two probability densities on X and Y , respectively. Then there exists 
a reflector which is a weak solution of the reflector problem. 

The proof is obtained in two steps. First the problem is solved in the case when 
the right hand side in ffTOj) is a finite sum of Dirac masses. In this case a constructive 
minimization procedure together with an apriori two-sided C° estimate of p is used 
to obtain the weak solution, which is also unique. The general problem is solved 
by approximating the right hand side in ( ITOi) by finite sums of Dirac masses and 
obtaining the solution i? as a limit of a sequence of special solutions Rk, k = 1,2, 
constructed on the previous step. The measures G{Rk, ■) are weakly continuous 
and converge to G{R, ■). Consequently, R is indeed a weak solution of the reflector 
problem. 

The procedure we have just described proves, in fact, that the reflector prob- 
lem admits a solution if the right hand side in fIlOp is any nonnegative probability 
measure on Y, possibly with a singular part. 

Existence of regular solutions was studied by X.-J. Wang [21] and P.-F. Guan 
and X.-J. Wang [lO]. 
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In the framework of Monge-Kantorovich theory the reflector problem was studied 
by T. Glimm and V. Oliker [8] and X.-J. Wang [22]. The following result was proved 
in [8], Theorem 4.1. 

Theorem 2.2. Let R he a weak solution of the reflector problem with the reflector 
map ^. Then ^^fi = v (that is, C, pushes fi forward to u) and it is a minimizer of 
the problem 



Furthermore, any other minimizer of fTS\) is equal to ^ almost everywhere on the 
set {x G I /(x) ^ 0}. 

Remark 2.3. Theorems \2.1\ and \2.Sl together, imply existence of minimizers to the 
problem fT^) . In addition, if I > on X and X is connected, these results imply 
that, except for a set of measure zero, any minimizer of ( flgj) (with fi and v as in 
ni]) ) is a reflector map associated with a closed convex reflector in R''"'""'^. Such 
reflector is unique up to a constant multiple of the function p(x) in see 

On the other hand, the minimization problem f|T3l) is a variant of the Monge 
problem on S*^ (see the beginning of section [3|, below). By a different method, 
in the framework of Monge-Kantorovich theory, the existence and uniqueness of 
minimizers to f|T3|) was proved in [8j and [22j under the additional condition that 
spt(/i) n spt(z/) = 0. Furthermore, if p and p are the functions in then Mo(x) = 
logp(x) and fo(y) = — logp(y) maximize 



over the set of pairs (u, f ) G C(X x Y) satisfying m(x) + v{y) < — log(l — x ■ y) for 
all (x, y) G X X Y. 

3 Background on the Monge-Kantorovich theory 

As observed in section [2], the reflector problem can be stated as a variant of the 
problem of Monge with the cost function c(x, y) = — log(l — x ■ y), x, y G S*^. 
In this section we recall some facts of the Monge-Kantorovich theory that will be 
needed to study an analogue of the reflector problem with general cost functions. 

We first fix some notation. For a set Z C R'^^"'^ we denote by V{Z) the set 
of Borel probability measures on Z. As usual, if Z is a closed subset of R'^"^^ and 



inf { / - log(l - X ■ ^'(x))rf/i(x) I = u} 



(13) 



x 
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7 e V{Z) then the support of 7 is the smallest closed set spt7 C Z such that 

7[ spt7] = l[Z] = 1. 

Suppose, X, Y C R'^^"'^ are closed sets. If /i G 'P(X), G ^(Y), we denote 
by r(/i, z/) the set of joint measures 7 on R'^"'"^ x R*^"^^ that have /j, and u as their 
marginals: ii[U] = -f[U x R'^+i] and 7[R'^+i x [/] = for Borel U C R'^+^ In 
fact, if 7 e r{ii, u) then spt7 C X x Y. 

Assume that we are given two probability measures fi and u on R^+^. Let T{fi, v) 
be the set of joint measures 7 on R'^^^ x R'^+^ that have and v as their marginals. 
Kantorovich's problem is to minimize the transport cost 

C(7) := j c(x,y)d7(x,y) (14) 

for some given c among joint measures 7 in r(//, i/), to obtain 

inf C[7]. (15) 
7er(/i,i/) 

Let T(/^, I/) be the set of Borel maps T : 'R!^^^ R''^^ that push ji forward 
to V : fj.[T-^{B)] = u[B] for all Borel sets B C R'^+^ The Monge problem is to 
minimize 

X[T] = [ c(x,Tx)d/^(x) 

over the set T(//, u). 

There is a natural embedding which associates to T G T{fi.u) a 77- := (id x 
T)#/x G r(^, z/), where id : R^+^ ^ R''^^ is the identity map. Since C[7t] = I[T] 
we conclude that 

inf Ch] < inf jm. (16) 

Throughout this section we use the notation R U {+00} = R and assume that 
c : R'^^^ X R''^^ R. We endow R with the usual topology so that c G C(R''+^ x 
R'^+^R) means that 

lini c(x,y) = c(x,y). 

(x,y)^(x,y) 

In particular, if c(x, y) = +00 then c(x, y) tends to +00 as (x, y) tends to (x, y). 

Definition 3.1. A subset S C R'^^^ x R''^^ is said to be c-cyclically monotone if 
for every natural number n and every {(xj,yj)}"^^ d S we have 

n n 

Xlc(xi,yi) < ^c(xi,y^(i)), 

i=l i=l 

for all permutation a of n letters. 
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This notion of c-cyclical monotonicity was introduced by Knott and Smith |13j . 
When c(x, y) = |x — yp, c-cychcal monotonicity is simply called cyclical mono- 
tonicity [ISJ. 

Proposition 3.2. Assume that X,Y C R'^"'"^ are closed sets, that ^ G V(X.), 
V G 'P(Y) and that c > zs lower semicontinuous on X x Y. Then, 
(i) there is at least one optimal measure 7o G T{fi, u). 

(a) Suppose that in addition c G C(X x Y,R). Unless C = +oo throughout 
r(/i, u), there is a c-cyclically monotone set S C R'^^"'^ x R"^"*"^ containing the support 
of all optimal measures in T{fi, u). 

Part (i) of proposition 13.21 can be found in [11] (Theorem 2.19). Note that that 
proof is interesting only in the case C ^ +oo. Part (ii) was established in [1] and [6] 
for c G C(X X Y, R). One can readily adapt the proof of Theorem 2.3 and corollary 
2.4 of [6] to cost functions c G C(X x Y, R). 

Definition 3.3. Suppose that ip : R'^+-^ — >■ R U {— C)o} is not identically — oo. Then 

(i) ip is said to be c-concave if there exists a set A C R'^^"'^ x R such that 

V;(x)= inf c(x,y) + A, (x G X). (17) 
{y,A)e.A 

(ii) The c-superdifferential d'^ip of ip consists of the pairs (x, y) G R'^^^ x R'^^"'^ for 
which 

c(x,y)-^(x) <c(v,y)-^(v), (18) 

for all V G R'^+^ 

(Hi) We define d'^xjj^K) C R'^+^ to be the set of y such that (x, y) G d'^tp^K). If 

E C R'^^^ d^iplE) is the union of the d''ilj{x) such that E. 

(iv) The c-transform of ip is the function y — >• ip'^iy) = infxgx{c(x, y) — ip(x.)}. 

Remark 3.4. Suppose that ip : R'^^^ ^ RU {— C)o} is not identically — oo and is 
given by [T7\). We have 



(i) ip^{y) > —A > — oo if{y, A) G ^ where A is the set in p7| ). Hence, ip^ ^ — oo. 

(ii) ip'^'^ = ip. 

Proof: The proofs of these remarks are well documented when c : R'^^^ xR'^+^ R. 
We verify that the same proofs apply when c may take the value +oo. 

If (y, A) G ^ then —A < c(x, y) — ip^x.) for all x G X. Hence ip'^{y), the infimum 
of c(x, y) — '?/'(x) over X, is not smaller than —A. This proves (i). 



11 



The inequality ifj'^ > ip which holds for general functions is readily checked. 
It remains to prove that when ( JT71) holds, then ip^'^ < ip. Fix x G X and let 
{(yn? ^n)}'^=i C ^ be such that 

ip{x) = lim c(x,y„) + A„. 

ra— >+oo 

By (i), A„ > -^'''(yra) and so, 

ijj{x.) > lim sup c(x, y„) - ^''(y^) > ^''''(x). 

ra— >+oo 

This proves (ii). QED. 



4 Existence and uniqueness of optimal maps 

Throughout this section, we assume that 
(A 1) ce Ci(R^+i X R^+i \ A). 

(A 2) c : R'^+^ X R'^"*'^ — > [0, +oo] is lower semicontinuous. 
(A 3) for any Xq G R'^^"'^, c(xo,Xo) = +oo. 

We are interested in probability measures /x, v for which C ^ +oo throughout 
r(/i,i/). Assume that c G C(R'^+^ x R'^+i \ A,R), where A := {(x,x) | x G R'^+^j 
denotes the diagonal. Proposition 14.11 provides a sufficient condition which ensures 
that C ^ +CX) throughout r(/i, u). Before stating that proposition, let us introduce 
the sets 

S'(a, 6) = {x = (xi, ■ ■ ■ , Xrf, Xd+i) I a < Xd+i < 
for -1 < a < 6 < 1. 

Proposition 4.1. Suppose that /i, z/ G ^(S'^) are Borel measures which vanish on 
{d — 1) -rectifiable sets. Then there exists 7 G T{fi, v) and e > such that 

|x-y| > e 

for all (x, y) G sj»i7. 

Proof: Since /z and vanish on {d — l)-rectifiable sets, the functions 

t ^ /i[S'(a, t)], t ^ u[S{a,t)] are continuous. (19) 
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Case 1. Assume first that there exists c G (—1, 1) such that 

spt/i C S'(-l,c), sptz/ C ^(c, 1). (20) 
Thanks to (I19p we may choose ei, £2 > such that 

-1 + ei)] = u[S{c, c + 62)] = 1/2. 

By (EOD -1 + ei < c. Set 

7 = 2(/i~ ® + /i+ ® 

where 

= ^|5(-l -l+ei)> /^"^ = /^|5(-l+ei,c), ^ = '^|s'(c,c+e2)5 '^'^ = i^|5(c+e2,l)- 

Note that 7 G r(/i, z/) and 

|x — y I > min{c + 1 — ei, 62} > 
for all (x, y) G spt7. This proves the proposition in this special case. 

Case 2. Assume that sptyU and sptz/ are arbitrary. We use f|T9l) to choose 
c G (—1, 1) such that 

fi[Sic,l)] = u[Si-l,c)]:=m. (21) 

If m = then 

sptyU C 5'(— 1, c), sptz/ C S{c, 1) 

and so, we reduce the discussion to the case 1. Similarly, if m = 1 we reduce the 
discussion to the case 1. 

Assume in the sequel that < m < 1. Set 

By (I2T1) . ^ and ^ are probability measures. They satisfy 

sptz/^ C S'(— l,c), sptyU^ C 5(0,1). 

Having that z/", //+ satisfy the assumptions of case 1, we may find 7 G r(^, ^) 
and e > such that 

|x-y|>e (22) 
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for all (x, y) G spt7. Similarly, there exists 7 G T^j^, jz:^) and e > such that 

|x-y|>e (23) 

for all (x, y) G spt7. 
Set 

7 = m7 + (1 — m)7. 
Then 7 G r(/i, i^) and by ([22]) -([23]), we have that 

|x — y| > min{e, e} > 

for all (x, y) G spt7. This proves the proposition. QED. 

Lemma 4.2. Suppose that c satisfies (A 1) - (A 3) and that X, Y C R'^^^ are 
closed sets. Suppose that S* C X x Y zs c-cyclically monotone and contains two 
pairs (xo,yo), (xo,yo) such that Xq 7^ Xo and jo 7^ fo- Then, there exists a function 
F G C(X X Y) depending only on the pairs such that 

) + c(x,y)<F(x,y) (24) 

for all (x, y) G S. 
Proof: Define 

F(x, y) = min{c(x, y^) + i?i(y), c(x, y^) + ^2(y)} 
Ri{y) = min{c(xo,y) + c(xo,yo),c(xo,yo) +c(xo,y)} (25) 
/?2(y) = min{c(xo, y) + c(xo, y^), c(xo, y^) + c(xo, y)} 

We use that Xq 7^ Xq and (A 1) - (A 3) to obtain that Ri,R2 G C(Y). This, 
together with the fact that jo 7^ fo yields that F is continuous on X x Y. If (x, y) 
is another element of S then setting 

(xi,yi) = (xo,yo), (X2,y2) = (xo,yo), (x3,y3) = (x,y) 
and using the c-cyclical monotonicity of S, we obtain f[24]) . QED. 

It is well known that a set is cyclically monotone if and only if it is contained in 
the subdifferential of a convex function [18]. An analogue of this result was proved 
by Smith and Knott [20j for general cost functions c : X x Y — > R. The following 
Lemma is a further extension that is needed to deal with cost functions satisfying 
(A 1) - (A 3) (and may be = +00 somewhere). Below, we check first that the proof 
in [19] extends to such cost functions and then we show that the infimum in f[26l) 
can be performed over the subset of y such that |x — y| > 5 > 0. 
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Lemma 4.3. Suppose that c satisfies (A 1) - (A 3) andX, Y C R'''^^ are compact 
sets. Suppose that S* C X x Y is c-cyclically monotone and contains two pairs 
(xcYo), (xo,yo) such that {xo,yo} Then, 

(i) S is contained in the c-superdifferential of a c-concave function ip : R^+^ — »• 
R U {— C)o} such that there exists 6 > satisfying 

V'(x) = inf {c(x, y) - ^'^(y) | |x - y| > ^j, (x G X). (26) 



(a) If (x, y) G d'^i/j then ip^x.) > — oo and |x — y| > 6 for some 6 > that depends 
only on ip. 

Proof: The expression ip{^) = infygY|c(x, y) — '?/^'^(y)|, x G X, is well-known in 

the literature. The only new and useful fact we want to point out is that V'(x) will 
be obtained by minimizing c(x, y) — ip'^iy) not on Y (as it is usually done), but on 

G Y I |x — y| > 6}. For completeness, we give the detailed proof below. 

Since X x Y is compact, the function F defined in 0241) attains its maximum. 
We use Lemma 14.21 and the fact that c is lower semicontinuous and equals +oo on 
A to conclude the following: if S" C X x Y is c-cyclically monotone and contains 
(xo,yo), (xo,yo) then there exists 6 > such that |x — y| > 26 for every (x, y) G 5". 
In particular, x 7^ y for (x, y) G S'. In particular, there exists ^5 > such that 
|x — y I > 26s for every (x, y) G S. 

As in [19j, we define 

{n— 1 n s 

c(x, y„) + c(xj+i, Yj) - c(xj, yj) \. (27) 
j=0 j=0 ^ 

Since c is finite on S and nonnegative on R'^"''^ x R'^"'"^, we conclude that ^/'(x) is 
well-defined. 

The c-cyclical monotonicity of S gives that 

n—l n 

c(xo, y„) + Y c(xj+i, Yj) - Y Yj) > 

and so, '?/'(xo) > 0. Taking n = 1, xi = Xq and yi = yo in (!27|) gives that V'(xo) < 0. 
We conclude that '?/'(xo) = and so, ip is not identically —00. This, together with 
the fact that ifj is clearly of the form (fTTl) yields that ifj is c-concave. 
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Claim 1. We have that 5" C d'^ip. 

Fix (x, y) e 5* and for each m integer, let {{^i,'yi)}"=i C S he such that 

lim 1pm = ^(x), 

m— >+oo 

where 

V'm := c(x,y„^) + c(xj+i,yj) - ^c(xj,yj). 

i=o i=o 

Setting (x„,„+i,y„„+i) = (x,y) we have that 

Tim Tlm.~\~^ 

-c(x,y)+V^„ = c(v,y„„+i)-c(v,y)+^c(xj+i,yj)- ^ c{xj,yj) > ^(v)-c(v,y). 

j=0 j=0 

Letting m go to +cx) we conclude that 

7/>(x)-c(x,y) >^(v)-c(v,y), (28) 

which proves claim 1. 

Claim 2. Whenever (x, y) G d'^ip, we have that ^/'(x) > — cxd and x 7^ y. 

Recall that (x, y) G d'^ip is equivalent to (|28l) . Setting (x, y) = (xo,yo), v = Xq 
in fl28l) . using the facts that '?/'(xo) = and Xq 7^ yo we obtain that ^^(xo) is finite. 
Next, if (uo, Vo) G 5, setting v = Uq we have that 

c(x, y) - ^(x) < c(uo, y) - V^(uo), (29) 

If y 7^ Xq we set = x^ in fl29|) to obtain the claim. If y = x^, we set = Xq and 
we use the fact that ipl'Xo) is finite to obtain the claim. 

Claim 3. The set d^ip is c-cyclically monotone. 

For the sake of completeness, we reprove this claim although it is a repetition of 
a known argument in [19]. If {(xj,yi)}^^i C d^ip, setting (x„+i,y.„+i) = (xi,yi) we 
have that 

c(xi, Yi) - ip{xi) < c(xi+i, yi) - ^/'(xi+i). (30) 
By claim 2, each term in fl30|l is finite and so, 

n n 

= ^^(xi+i) -^(xj) < ^c(xi+i,yi) -c(xi,yi), 

i=l i=l 
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which proves the claim. 

We use Lemma |472| the facts that d'^%lj is c-cychcally monotone, that (xo,yo), 
(xo,yo) £ d^'ip, that Xq 7^ Xq and jo 7^ Yo, to obtain the existence of some 5 > 
such that |x — y| > 25 for all (x, y) G d'^ip. By remark [33] (ii) -i/^ = [ip'^Y and so if 
X G X, there exists a sequence {ynjj^i C Y such that 

7/>(x)= lim c(x,y„)-^^(yn). (31) 

Since Y is compact, we may extract from {ynj^JLi a subsequence (which we still 
label {ynjJ^Li) that converges to some y G Y. Recall that the function ip'^ is upper 
semi continuous as an infimum of upper semicontinuous functions and therefore, fl31l) 
yields 

7/'(x) > c(x,y) -^/'^(y). 

This proves that V^(x) = c(x, y) — ip'^iy) and so, (x, y) G d^ip. Hence |x — y| > 25 
and y is a minimizer in (|26|1 . QED. 



Remark 4.4. Suppose that fi G V(X.) and v G 'P(Y) /iawe no atoms and that 7 
minimizes C over r(/i, z/) anc? t/iai ^(7) < +00. Then, 7[A] = and so, spt{'j) \ A 
contains at least one element, say (xo,yo). Also, 'y[E] = where 



E = ({x„,yj X y) U (^X X {x„yj) 



Hence, the set X x Y\ (i?U A)) is nonempty, and so, it contains an element (xo,yo). 
Note that x^, y^ ^ {x^, yo}. 

We now further specialize the set of cost functions under consideration by as- 
suming that 

( ^ /^(^) if x^y 
c(x,y) = < , (32) 
I +00 if X = y 

where / G C^(0, +00) is such that 

limt_o+ l{t) = +00 

r(t)| >o (t>o) ^ ^ 

Define 

git)=ti2-t)int)r (34) 
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Note that if g is monotone on (0,2] then 5'[|,2] is a closed interval on which g ^ 
exists and is continuous. Define 

M(a,x)= [l-g-\\8i\') 

and the closed set 



Note that M is continuous on Ks for all 5 > 0. 

When X G S'^ we next denote by VgdC the tangential gradient of x — > c(x, y) at 
X G S'^ and let Tx be the tangent space to at x G S^. 

Lemma 4.5. Assume that c is given by [3^) such that I satisfies ^33\) . Assume 
that g is monotone on (0,2]. If ^ a. & R'^'^-'^, x, y G S*^ and VgdC(x, y) = a then 
y = M(a,x). 

Proof: Recall that if x G S'^ then the orthogonal projection of y G R^+^ onto Tx is 
y|[ = y — (x ■ y)x. Hence if y G S** then 

|y|l|2 + (x-y)2 = l (35) 



Setting a := VgdC(x, y) yields that 



-^'(^V^)yii- (36) 



2 

2 ■ 



This, together with (135|) and the fact that 1 — x • y = for x, y G S*^, yields 

that 

Thus, 

l-x-y = ^-i(|ap). (37) 
We use f p6l) . fl37|) and the fact that 1 — x ■ y = ^^^^^ to conclude that 

y = (x-y)x + y|| = fl-g-\\si\^))^-^^ ^ = M(a,x). 

QED. 
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Theorem 4.6. Assume that c is given by ^3B\) with I G C^(0,+oo) and satisfying 
Assume that g given in [34^ is monotone. Let yU, G 7^(8*^) be Borel measures 
on the sphere S*^ that vanish on {d — 1) -rectifiable sets. Then 

(i) there exists a unique measure 7o that minimizes C over T{fi, u). 

(a) There exists a unique map Tg : S*^ ^ that minimizes I over T{^.,v). 
Furthermore the essential infimum infx |ToX — x| > 0. The is uniquely determined 
and coincides with the measure (id x To)#fi. 

(Hi) The map To is invertible except on a set whose fi measure is null. 

Proof: (i): Existence of tlie measure jo- By proposition HTTl tliere exists 7 G r(/i, u) 
and e > sucli tliat 

|x-y|>e (38) 

for all (x, y) G spt7. Since |x — y| < 2 on S*^, ( l38i) and the fact that c satisfies 
(A 1) ensures that c is uniformly bounded on spt7. This proves that C['j] < +00. 
By proposition 13.21 there is at least one optimal measure 70 G T{fj,,u) and there is 
a c-cyclically monotone set S C R'^^^ x R'^+^ containing the support of all optimal 
measures in r(/x, z/). By remark the support of 7o and hence S contains pairs 
(xo,yo), (xo,yo) such that Xq, yo {xo,yo}. Lemma H73] ensures existence of 5 > 
and a c-concave function such that S is contained in the c-superdifferential of 
and 

^(x) = inf|c(x,y)-^^(y)| |x-y|>5}, (x G X). (39) 

Clearly i/j is continuous and so is ip'^. Hence, d'^il){^),d'^ip'^{y) 7^ for all x, y G S*^ 
and 

(x,y)Ga^^ ^ V^(x)+^^(y)=c(x,y) (40) 

We use again fl39l) and the assumption that / G C^(0, 00) to conclude that ip is 
semi concave on S*^ and so, it is differentiable everywhere except for a set iV C S*^ 
which is {d — l)-rectifiable. 

Proof of (ii): existence and uniqueness of To] 7o = (id x To)#^; uniqueness of 70. Let 
X G S'^ \ and y G d'^ip{'x). We use the fact that the function t —>■ ip{t) + ip^{y) — 
c(t,y) is difi^erentiable and attains its maximum at x to conclude that 

Vs<iV^(x) = V^.c(x,y). (41) 

This, together with the Lemma [4.5[ implies that y = M(Vsd^/'(x), x). Since 

|x-y|>5, (42) 
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we obtain that (Vsd'?/'(x), x) G K^. Hence the map To defined by 

T„x:=M(Vsd^(x),x) (43) 
is a Borel map. Note that by fH2]) . the map To satisfies 

|ToX — x| > 5 

for all X e S'* \ A^. 

Because /i vanishes on N we have that 7o[A^ x S*^] = and so, To is defined ix 
almost everywhere. Note that we have proved that 

d^'^l) \{Nx S*^) C graphTo. (44) 

Since 7o[A^ x S'^] = and 7o G r(/i, u), we obtain that 

j L(x, y)rf7(x, y) = j Tox)d7o(x, y) = L(x, Tox)dfi{x), (45) 

for all L G C(R'^+i x R'^+i). By (jlS]) we have that 

7o = (id X To)#/i, T#/i = i^. 

This proves that 7o is uniquely determined. We use f l45l) with c = T to obtain that 
C[7o] = X[To]. By f ITBl) we conclude that To minimizes T over T(yU, z/). 

Furthermore, if Ti is another minimizer of X over T(/i, z/) then 7o = (id x Ti)#/i 
and (x, Tix) G spt7o C d'^ip for almost every x. By f H4|) we have that Ti (x) = ToX 
for /i almost every x. This proves that To is uniquely determined. 

Proof of (iii): invertibility of the map Tq. The analogue of (1391) for -ip^ gives that ■ip'^ 

is semiconcave and so, the set N where ip'^ is not differentiable is {d— l)-rectifiable. 
Substituting ^ hj u , the above reasoning yields that the map 

5oy = M(Vsd7/;^(y),y), (y GS'^\iV), 

is such that So#i^ = fi, and 

d'^tfj' \{Nx S*^) C graph^o- (46) 

We use dS]) , (UnD and the fact that d^ip = d^ip" to conclude that id = 5*0 o To on 
S*^ \ (A^ U S~^{N). Since is (c? — l)-rectifiable and S'o#z/ = yU we have that 

4S;\N)] = u[N] = 0. 

Thus, To is invertible on S'^ up to a set whose /i measure is null. QED. 
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Remark 4.7. If it is assumed only that I G C^(0, oo) and fi is absolutely continuous 
with respect to the standard measure on S** then it is easy to see that ip, ip'^ are 
differentiable a.e. on S*^. Indeed, it follows from ^3^, ^'^^ the assumption 
that I G C^(0, cxo) that ip,ip^ are locally Lipschitz. Then by Rademacher's theorem 
ip,ip'^ are differentiable a.e. Note, however, that this result is weaker than the one 
established in Theorem \4.6[ 
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